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ABSTRACT 

Aims. In the Local Group spiral galaxy M 33, we investigate the correlation between the star formation rate (SFR) surface density, 
S SFR , and the gas density S g!ls (molecular, atomic, and total). We also explore whether there are other physical quantities, such as the 
hydrostatic pressure and dust optical depth, which establish a good correlation with S SFR . 

Methods. We use the Ha, far-ultraviolet (FUV), and bolometric emission maps to infer the SFR locally at different spatial scales, and 
in radial bins using azimuthally averaged values. Most of the local analysis is done using the highest spatial resolution allowed by gas 
surveys, 180 pc. The Kennicutt-Schmidt (KS) law, Ssfr x XL, is analyzed by three statistical methods. 

Results. At all spatial scales, with Ha emission as a SFR tracer, the KS indices are always steeper than those derived with the FUV 
and bolometric emissions. We attribute this to the lack of Ha emission in low luminosity regions where most stars form in small 
clusters with an incomplete initial mass function at their high mass end. For azimuthally averaged values the depletion timescale for 
the molecular gas is constant, and the KS index is » H , = 1.1 ±0.1. Locally, at a spatial resolution of 180 pc, the correlation between 
Esfr and E gas is generally poor, even though it is tighter with the molecular and total gas than with the atomic gas alone. Considering 
only positions where the CO 7=1-0 line is above the 2-cr detection threshold and taking into account uncertainties in E H , and S SFR , 
we obtain a steeper KS index than obtained with radial averages: «h 2 = 2.22 ± 0.07 (for FUV and bolometric SFR tracers), flatter than 
that relative to the total gas (ttH 10t - 2.59 ± 0.05). The gas depletion timescale is therefore larger in regions of lower Esfr- Lower KS 
indices («h 2 = 1-46 ± 0.34 and «h 2 = 1.12) are found using different fitting techniques, which do not account for individual position 
uncertainties. At coarser spatial resolutions these indices get slightly steeper, and the correlation improves. We find an almost linear 
relation and a better correlation coefficient between the local Esfr and the ISM hydrostatic pressure or the gas volume density. This 
suggests that the stellar disk, gravitationally dominant with respect to the gaseous disk in M 33, has a non-marginal role in driving 
the SFR. However, the tight local correlation that exists between the dust optical depth and the SFR sheds light on the alternative 
hypothesis that the dust column density is a good tracer of the gas that is prone to star formation. 

Key words. Galaxies: individual (M 33) - Galaxies: ISM - Galaxies: Local Group - Galaxies: spiral 



■ 1. Introduction good correlation for the atomic and total gas densities (KS in- 

' dex «h, m = 1 .3 + 0.3), but not fo r the molecular gas density. One 

J> . The gas-to-star conversion process is one of the most important decade later. iKennicuttl i 1 998bl) studied the relation between the 



ingredients for galaxy evolution. The rate at which stars form total gas surface density and the Ha SFR density averaged over 
in a galaxy at a given epoch depends not only on the available the entire galaxy, using observations for 61 normal spiral and 36 
gas reservoir but also on the abili ty of the gas to colla pse and starburst galaxies. Due to the rather large scatter in the KS re- 
fragment. The first seminal papers dSchmidtll 195911 1963b related lation, the KS index was highly dependent on the method used 
the star formation rate (SFR) to the atomic gas densities using a to fit the data. For the 61 normal galaxies, a least-squares fit on 
power law. Comparing the gas density w ith the number o f young the SFR density yields n H|ot = 1.29 + 0.18, while a bivariate 
stellar objects in the solar neighborhood Schmidt] dl959h derived least-squares regression leads to n HlM = 2.47 + 0.39. Likewise, 
a power law index n = 2, i.e. Esfr oc Ej^ s . Va lues of n ap- these two methods to estimate the KS indices for the sample of 
proxi mately 1.5 to 2.0 were further confirmed bv lGuibert et alj 36 starburst galaxies lead to values of n H , , = 1.28 ± 0.08 and 
(119781) . using more precise data on the radial and vertical dis- n Htm = 1.40 + 0.13, respectively. The better agreement between 
tributions of the interstellar gas and a variety of young stellar the two fitting methods mainly reflects the higher gas and SFR 
objects. dynamic ranges (three orders of magnitude) spanned by the sam- 
Papers on external galaxies often use the term Schmidt law P le of starburst galaxies with respect to the sample of normal 
or Kennicutt-Schmidt (KS) law to relate the surface density of galaxies (one order of magnitude). The derivation of the KS in- 
gas to the SFR per unit surface area, since these quantities are dex 18 ver y sensitive to the method used to fit the data when there 
the observables for external galaxies. For half a ce ntury, nu- ls not a wlde dynamical range of the variables. 
merou s studies have been done on the KS l aw (s ee IKennicuttl Today several questions regarding the KS law are still open. 
1998a, for areview). For instance. iKennicuttl d 1 9891) studied how First, does the SFR surface density correlate better with the total 
the globally averaged SFR in a galaxy traced by Ha emission (atomic plus molecular) gas surface density, as suggested by the 
correlates with the mean atomic, molecular, and total (H 1+H2) globally averaged studies, or only with molecular surface den- 
gas surface densities. For a sample of 15 galaxies they found a sity, since stars condense out of molecular gas? Second, which 
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is the best tracer to characterize the SFR? Is Ha a good local cur- 
rent SFR tracer? Does the incompleteness of the IMF at the high 
mass end in regions w here only small cluste rs form make Ha an 
unreliable SFR tracer dCorbelli et al.ll2009h ? The far-UV radia- 
tion traces recent SFR averaged over a longer period of time (up 
to 100 Myr): how well does it correlate with the gas density lo- 
cally? On which scale and at which wavelength does the infrared 
radiation trace star formation? How much of the evolved stellar 
population is contributing to dust heating? 

Third, the dependence of the KS law on the spatial scale 
considered is not yet clear. It is of interest to study at which 
scale the KS law might break down and how the power law 
index n changes as we sample smaller and smaller regions. A 
first step towards a spatially resolved KS law was performed us- 
ing azimuthally averaged v alues of SFR and gas densities (e.g. 
Martin & Kennicuttl l200lt IWong & Blitzl 120021 iBoissier et all 
2003112007b . For instance. fWong & Blitz (2002) studied a sam- 



ple of molecule-rich spiral galaxies and found that the SFR den- 
sity is more strongly correlated with the H2 surface density than 
with the total gas surface density. They derived «h, = 1-2 with 
an average value of 1 .4 using the molecular gas surface density 
and a radially varying extinction correction. For the total gas sur- 
face density the average KS index found by the same authors is 
steeper, h Hm ~ 1.7. 

As the resolution of the telescopes increased, the KS law 
has been examined locally at severa l spatial resolutions (e.g. 
ICalzetti et aLll2005t iBigiel et alj|2008l) . However, feedback pro- 
cesses linked to SF activity cast doubt on the applicability of 
the KS law on very small scales. Photodissociation or photoion- 
ization radiation from massive stars or stellar winds can locally 
remove the molecular hydrogen and make H n regions not spa- 
tially coincident with peaks of CO emission. This effect is visible 
in M 33 for example when comparing t he Ha emission map of 
this galaxy with C O 7=1-0 line maps (Engarg iola et al.| [2003; 
lHeveretalJ2004h .M33 is a galaxy with a generally low extinc- 
tion, and the recent Spitzer maps of this galaxy at 8 and 24 pm 
emission co nfirm that the disp l acement is not due t o embedded 
H n regions (IVerlev et alj|2009t ICorbelli et al]|2009l) . 

In order to addr ess these questions, rec ent space mission data 
in the UV (GALEX , \Gil de Paz et al.l2007T) and in the IR (Spitzer, 
IWerner et al.l l2004) have been used together with ground-based 
maps in the Ha recombination line. Such a multiwavelength 
database enables accurate diagnostics of the SFR at high angu- 
lar resoluti on in nearby galax ies. The H 1 Nearby Galaxy Survey 
(THINGS. IWalter et al.ll2008l) has provided H 1 and CO maps for 
many nearby galaxies at an angular resolution of 500-800 pc, 
thus offering the possibility to examine the relation between the 
gas sur face density and S FR density locally, on smaller spatial 
scales. Bigiel et al. (2008) find that the best correlation is be- 
tween the H2 surface density and the SFR (traced by a combi- 
nation of the far-UV and IR surface brightness). The average 
KS index is: 1.0 + 0.2. The value «h, = 1 implies that the SFR 
scales with the mass surface density of molecular gas, i.e., that 
the timescale and the efficiency of SF (the fraction of gas mass 
converted into stars) is constant. 

The THINGS result of unit slope differs from what has been 
found in other nearby galaxies by other groups using multi- 
wavelength data. Thilke r et al.l d2007l) derive a KS index of 1 .64 
(1.87) using £h 2 (£h[+h 2 ) m molecular dominated regions of 
NGC7331. In M51 a KS index of 1.56 has been found by 
iKennicutt et al.l (120071) for the spatially resolved SFR density 
on 500 pc scales. The index refers to the total gas density and 
has been derived by taking into account only bright H n regions 
in the center or along the arms of M 5 1 . The KS index for the 



H2 surface density in M51 is slightly lower, 1.37, very similar 
to that found for azimuthally- averaged quantities in the BIMA 
samp le or previ ously in M 33 fWp ng & Blitz 20Q2 iHever et all 
120041) . In M 3 1 dBraun et alj|2009l) the best correlation at 1 1 3 pc 
scales is established between the SFR density and the total gas 
density, with a KS index similar to that found for M 5 1 . The dif- 
ferent power law indices found in the literature, 1 < n < 3.5, 
imply that either the SF law differs from galaxy to galaxy (or 
from region to region such as arms versus whole disk) or that 
the KS law index is very sensitive to the method used for deriv- 
ing it (spatial scale, SFR indicator, extinction corrections, CO to 
H2 conversion, background subtraction, etc.). 

Due to its proximity, large angular size, and rather low in- 
clination, the Local Group spiral galaxy M 33 is a unique tar- 
get to investigate the physics underlying the KS law in a late- 
type galaxy. M 33 has about 1/3 of its baryons in gase ous form, 
a small molecular fractio n and a low dust abundance (Corbelli 
120031: IVerlev et al .1 120091) . Early attempts to test the KS law in 
M 33 considered only the atomic gas because of its low molec- 
ular content, given the limited spatial resolution and sensitiv- 
ity of molecular gas observations. In the early 1970's, using 
stars and Hn r egion counts together with neutral hydrogen gas, 
[ Madore et ail (1 19741) . found a KS index « H , ~ 2.35 + 0.26. 
Newton ( 1980) noticed that the density of H 11 regions on scales 
of ~ 300 pc had a weaker dependence on the H 1 surface den- 
sity in the inner regions of the galaxy than in the outer ones. 
More recently, using the first unbiased census of 12 CO 7=1-0 
line emission in M33 and t h e far- infrared emission map pro- 
vided by IRAS, IHever et akl d2004l) found a strong correlation 
between the azimuthally averaged SFR density and the aver- 
age molecular gas surface computed for annular regions 250 pc 
wide. For the molecular gas surface density, the KS index is 
«h, = 1.36 ± 0.08, while for the total gas surface density, it 
is much steeper, tin m = 3.3 + 0.07. This steeper index is con- 
sistent with the molecular index, given the low molecular gas 
fraction in this galaxy which seems regulated by the balance be- 
tween the gas pressure (acting on the H2 formation rate) and the 
dissociation radiation dElmegreenlll993at IWong & Blitzl [2002; 
IHever et al.ll2004hlBlitz & Rosolowsky||2006l)~ 

Another feature of M 33 that makes it interesting for such 
a study is the inconsistency of the disk stability wit h the ongo- 
ing SF dKennicuttlll9"89t iMartin & Kennicuttl 1200 ll) . For many 
years most of the star-forming disk of M 33 has been known to 
be stable according to th e simple Toomre gravitational stability 
criterion (Toomre 1964) if only the gas su rface density is con- 
sidered (Elmegreen 1993a). Corbelli (2003) has shown however 
that the Toomre criterion predicts correctly the size of the unsta- 
ble star forming region of the M 33 disk when the stellar gravity 
is considered in addition to that of the gaseous disk. Apparently 
the stellar disk plays an important role in driving the disk insta- 
bilities which trigger SF. Its gravity compresses the gas and it 
can affect the SFR as well. If, for example, the density, or the 
free fall time, of a cloud depends on the disk gravity perpendic- 
ular to the plane, then the SFR density might not correlate with 
the gas density alone. Also if SF happens only in self gravitating 
clouds, then the presence of diffuse molecular material in CO all 
disk surveys might weaken the expected correlation. 

The wealth and quality of newly available data for M33 
has drastically increased during the pa st years thanks to re- 
cent space missions in the UV (GALEX, Gil de Paz et al.l l2007) 
and in the IR (Spitzer, IWerner et alj|2004b . The high resolution 
and sensitivity of the rec ent multiwavelength database for M 33 
dVerlev et al.ll2007l 12009) makes it an ideal target for investigat- 
ing the radial and local relations between the various SFR trac- 
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ers and g as densities. The spatial resolutio n of the CO 7=1-0 
survey dcorbeilill2003t iHever et al.H"2004 of M 33 is 45 arcsec, 
similar to the resolution achieved by Spitzer at 160 pm (~40 arc- 
sec). We can therefore test the KS law at a spatial resolution of 
180pc. To shed light on the physical basis of the scaling relation 
between SFR and gas density in late type spiral galaxies, we test 
whether a simple relation between SFR and gas surface density 
exists in M 33. We do this using several SFR tracers and spatial 
resolutions applying different fitting methods. Galaxies such as 
M 33, which do not have a high SFR per unit surface area and 
have a nearly constant gas surface density, might not establish 
a tight KS relation; hence the difference between various fitting 
methods can be large. 

This artic le is the fourth in a series d edicated to the SF in 
M33 after I Verley et all (120071) (hereafter IPaper J), IVerlev et all 
(120091) (h ereafter Paper III) and ICorbelli et al.1 (120091) (hereafter 
Pape r IlJ) . It is organized as follows: Section |2]presents the data 
and the methodology used to derive SFRs and gas surface den- 
sities. Section [3] compares the neutral gas distribution, atomic 
and molecular, with the Ha emission line map. In Sect. [4] we 
study the KS law considering azimuthal averages of different 
SFR and gas surface density tracers, and in Sect. [5] we examine 
the local KS law at various spatial resolutions, using two differ- 
ent methods. In Sect. [6] we investigate whether the SFR per unit 
area establishes a better correlation with the gas volume density, 
and we outline possible biases when using Ha emission as SFR 
tracer on a local scale and the CO line luminosities as molecular 
gas tracer. Our conclusions are summarized in Sect. [7] 



2. The data sets and methodology 

Here we describe the multiwavelength data set that has been 
compiled and also the photometric methods and SFR diagnos- 
tics. These diagnostics will be used both for the azimuthally- 
averaged KS law and for the local KS law on a series of increas- 
ing spatial scales. 

2. 1 . Ultraviolet and Ha line images 

To investigate the continuum ultraviolet (UV) emission of 

M33, we use Galaxy Evolution Explorer (GALEX) mission 

■ 1 1 1 

(Martin et al. 2005) data, in particular the data distributed by 
iGil de Paz et alJ d2007l) . A description of GALEX observations 
in the far-UV (FUV, 1350-1750 A) and near-UV (NUV, 1750- 
2750 A) relative to M 33 and of the data reduction and calibra- 
tion procedure can be found in lThilker et al.l d2005l) . 

To trace ionize d gas, we adopt th e narrow-line Ha image of 
M 33 obtained by Greenawalt d!998l) . The reduction process, us- 
ing standard IRAlfl procedur es to subtract the cont i nuum emis- 
sion, is described in detail in iHoopes & Walterbosl d2000l) . The 
total field of view of the image is 1 .75 x 1 .75 deg 2 (2048 x 2048 
pixels of 2'.'028 on a side). 

2.2. Infrared images 

Dust emission can be investigated through the mid- and FIR 
data of M 33 obtained with the Multiband Imaging Photometer 
for Spitzer (MIPS) instrument dWerner et al.1 120041 iRieke et al.1 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



2004])- The complete set of MI PS (24 70, and 160 pm) im- 
ages of M33 is described in IPaper It the Mopex software 
dMakovoz & Marleaull2005l) was used to gather and reduce the 
Basic Calibrated Data (BCD). We chose a common pixel size 
equal to 172 for all images. The images were background sub- 
tracted, as explained in lPaper j The spatial resolutions measured 
on the images are 6", 16", and 40" for MIPS 24, 70, and 160 pm, 
respectively. The complete field-of-view observed by Spitzer is 
very large and allows us to achieve high redundancy and a com- 
plete picture of the star forming disk of M 33, despite its rela- 
tively large extension on the sky. 



2.3. 21 -cm and millimeter data 

Several data sets are available to examine the atomic and 
molecular gas distributions: these include the Westerbork Radio 
Synth esis Telescope (WRST) array data dDeul & van der Hulstl 
Il987l 24" x 48" spatial resolutio n) and Arecibo single dish 
survey dCorbelli & Schneidel fl997l 4' sp. res.). For the molec- 
ular gas emission as traced by the CO 7=1-0 rotational 
line, we can use th e Berkeley Illinois Ma ryland Association 
(BIMA) array data dEngargiola et al.1 120031, 13" sp. res.) and 
the Five Colle ge Radio Astronomy Observatory (FCRAO) sin- 
gle dish data dCorbellil 120031: iHever et all 12004, 45" sp. res.) 
or the map obtained by com bining the two surveys as described 
by iRosolowsk v et all d2007l) . The interferometers in general re- 
cover less flux than single-dish data, since they tend to filter the 
diffuse emission from structures much more extended than the 
primary beam resolution. In fact, BIMA obser vations of M33 
conta in roughly half of the single-dish flux (Engargiola et al.l 
120031) : this implies that single-dish data are potentially more 
reliable to measure the total gas column density in M 33, even 
though it is not clear which one establishes a better correla- 
tion with the SFR. We shall use the WRST and the FCRAO 
data for the atomic and molecular gas distributions res pectively, 
which have a comparable spatial resolution. Following ICorbellil 
(2003), the CO measurements were converted to H2 mass col- 
umn densities (M pc~ 2 ), using the standard conversion factor, 
X = 2.8 x 10 20 cm" 2 (Kkm s -1 ) -1 . 



2.4. Star formation tracers 

We can use several SF tracers to test the KS law. Different trac- 
ers are sensitive to different timescales of SF, and their accu- 
racy to trace SF episodes strongly depends on the SF history and 
dust content of the galaxy under scrutiny. Ha emission traces gas 
ionized by massive stars in recent bursts of star formation over 
timescales of 10 Myr or so. The FUV luminosity corresponds to 
relatively young stellar populations (< 100 Myr), so can be con- 
sidered as complementary to Ha in terms of sensitivity to short 
timescales. 

To convert Ho- and FUV emission into SFR, we first cor- 
rect both H a and FUV for extinction using the formalism de- 
veloped by Calzettjl (12001 1) . This empirical approach relates 
the extinction to the TIR and FUV luminosities. The TIR flux 

in W m~ 2 pc~ 2 ) i s the total IR flux from 3 to 1000 pm 

Dale & Heloull2002l) defined as: 



F(TIR) = 10" 14 x [19.5 F v (24) + 3.3 F v (10) + 2.6 F v (160)] , (1) 
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where F„(24), F v ( 70), and F v (160) are the MIPS flux densities 
in Jy pc~ 2 (see alsoPagerllj). The extinction correction can then 
be written as: 



1 Z/TIR) 

A FUV = 2.5 x C x log I x — + 1 

s V 1.68 L(FUV) 



(2) 



The value of C is unity if considering star forming regions. 
Since we sample star forming regions and the I SM with older 
populations, we adopt an average value C — 0.7 (Pap er III) . For 
Ha extinction we shall use An a = 0.3 x Apuy Following the cal- 
ibrations given in iPaper Ifl, the SFR for the Ha and FUV wave- 
lengths are given by the following formulae: 



SFR(Ha) [M yr" 1 ] = 8.3 x 10~ 42 L(Ha) [erg s" 1 ] . 
SFR(FUV) [M G yr 1 ] = 8.8 x 10~ 44 L(FUV) [erg s" 1 ] 



(3) 
(4) 



Furthermore, Paper II has shown that in a galaxy of low 
dust content such as M 33, multifrequency diagnostics can help 
to determine the SFR more accurately. We adopt one "hybrid" 
SFR tracer in addition to Ha and FUV tracers: the bolometric 
luminosit y, approximated by a linear combination of the FUV 
and TIR (Thil ker et alj|2007l) . This is potentially a reliable SF 
tracer in regions which are more extended than individual Hn 
regions. It refle cts both t he unobscured and obscured regions of 
the galaxy (see IPaper lilt) and is sensitive to a variety of star- 
formation timescales. To estimate the b olometric luminosi ty, we 
used the following equation, given by Thilk er et al.l d2007l) : 



L(bol) = v FUV L v , obs (FUV) + (1 - 77)L(TIR) 



(5) 



We know that the interstellar radiation field which heats the 
dust has a contribution from the old stellar population. Some 
of the IR emission may also come from dust heated by evolved 
stars (AGB), a nd thus not be directly associated with recent SF 
episodes (see Paper II). This heating might dominat e for dust 
outside H n regions. The fa ctor rj accounts for this (Bell 2003; 
Igl esias-Paramo et al. 2006), and we shall u se r\ - 0.3 . The bolo- 
metric SFR was calculated according to lThilker et al.l d2007l) : 



SFR(bol) [Moyr 1 ] 



4.6 x 10" 44 L(bol) [erg s" 1 ] 



(6) 



This equation is bas ed on the calibration of 
Id esias-Para mo et all (|2006) derived from Starburst 99 models 
dLeitherer et alJ 119991) . with a Salpeter IMF (0.1 - 100 M Q ), 
solar metallicity, and a continuous SF. 

2.5. Radial profiles and aperture photometry 

In Sect. |4] we shall investigate the azimuthally averaged KS 
law by correlating the azimuthall y averaged values o f the SFR 
to gas surface density. Following iHever et al.1 ([2004), we com- 
pute Ssfr and 2 gas as the mean values within elliptical an- 
nuli spaced by 0.24 kpc and centered on the galaxy center 
(01 h 33 m 50 s .90, +30°39'35'/8). The annuli are assumed to be cir- 
cular rings viewed at an inclination of 54° and with the line of 
nodes at a position angle of 22.5°, thus representing the spa- 
tial or ientation of the M33 dis k with respect to the line of 
sight dMcConnachie et all 120061) . We performed this analysis 
on molecular, atomic and total gas tracers (from FCRAO and 
WRST data) and for each of the SFR tracers as described in the 
previous section. These azimuthal averages give £sfr and S gas as 
a function of radius, the central radius of each ring. We then cor- 
rect E gas and the radiation emitted from the newly formed stars 
for the disk inclination to obtain face-on values. 



In addition to analyzing azimuthal averages, we shall inves- 
tigate the local KS law using the highest spatial resolution pos- 
sible for our dataset (see Sect. [5J. This is the spatial resolution 
of FCRAO CO 7=1-0 map, similar however to that of 21 -cm 
and 160 /mi maps. In order to compare Esfr and E gas locally we 
first degraded the 24, 70, 160 yum, FUV, and Ha images to the 
lowest resolution, i.e. to 45", corresponding to the spatial reso- 
lution of FCRAO molecular map. This corresponds to sizes of 
180pc at the distance of M33. This is still large enough not to 
resolve out single molecular clouds: molecular cloud s sizes in 
M 33 are smaller than 100 pc dRosolowskv et al.ll2003b . Then, at 
every position observed in the FCRAO map, we performed aper- 
ture photometry on each of the images with a beam of 45" (full 
width half maximum). This gives a total of 7664 positions over 
the disk of M33. Finally, we converted the photometric results 
to the appropriate surface density units to examine the local KS 
law. 

To investigate how the KS law varies as a function of the 
spatial scale we evaluate the KS law on coarser spatial scales by 
averaging the surface densities over adjacent positions. 

3. Comparison of neutral and ionized gas 
distributions 

Figure[T]shows the H i contours (at logarithm levels of 0.55, 0.85, 
1.03, 1.15 M pc~ 2 ) overplotted on the Ha image. The Hi atomic 
gas distribution in M 33 is very filamentary, and the H n regions 
lie on the high surface density Hi filaments. But the bright Ha 
emission knots in M 33 are generally not coincident with the lo- 
cation of the neutral gas peaks. The Hi contours follow the Ha 
emission on larg e spatial scales, but n ot locally, on small scales. 
This agrees with Wright et al] d 1972b . who found no correlation 
between the positions of the ten brightest H n regions and the 
H i density peaks. In the northern spiral arm, south of NGC 604, 
there is relatively strong H i emission, but little Ha em ission as- 
sociated with it, as already noted by iNewtonl (Il980l) . The Hi 
emission extends very far from the center, much beyond the op- 
tical disk, with a slower radial decline than that of the Ha surface 
brightness. The H i contours even trace the northern plume where 
the Ha emission is faint. The maxima of the H i emission are in 
the southern half of the galaxy, exactly on the southern arm, and 
not in the geometrical center of M 33, where the gas mass den- 
sity is dominated by molecular gas, and several areas appear to 
be totally devoid of H i. 

We use the FCRAO CO 7=1-0 emission map to trace 
the H2 molecular gas, with the conversion factor as given in 
Sect. 12.31 In Fig. [2] the levels of the contours are, in logarithm, 
0.90, 1.11, 1.25, 1.36 M pc~ 2 . Above the detection threshold 
of FCRAO measurements (3<x ~ 0.8 K km s~'), molecular gas 
is visible only in the inner disk of the galaxy, up to about 3- 
4 kpc. CO contours closely follow the southern spiral arm of the 
galaxy. In the north is a displacement of the CO emission peaks 
relative to the arm as traced by H a. This displacement is also 
evident in the BIMA data of M 33 ( Engargi ola et al.ll2003l) . For 
the GMCs of BIMA. lEngargiola et all (12003) showed that there 
is clustering of GMCs and Hn regions out to a separation of 
150pc. But only 67% of the detected GMCs have their centroid 
position within 50 pc of an Hn region. Along the minor axis, 
on the western side of the galaxy, CO is absent, while there are 
bright H 11 regions visible in Ha emission as well as in the IR. 
No CO peaks are visible close to IC 133, although CO emission 
is detected around the bright H 11 region NGC 604. 

Unlike luminous spiral gal axies, the gas com ponent of M 33 
is dominated by atomic gas (IHever et al.l |2004|) . The atomic- 
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Fig. 1. Equal density H i contours (in logarithm: 0.55, 0.85, 1.03, 
1.15 M pc~ 2 ) superimposed over the Ha image. 
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Fig. 2. Equal density CO contours (in logarithm: 0.90, 1.11, 
1.25, 1.36 M pc~ 2 ) superimposed over the Ha image. 



to-molecular fraction is extremely high ( 1 to 10) compared t o 
molecule-rich galaxies (e.g., 0.05 to 0.4, [Wong & Blitz 2002). 
Small molecular c louds have been found beyond 3-4 kpc 
dGardan et al.l l2007). The paucity of molecular gas at large radii 
compared to the atomic gas, which has no radial decline, can be 
explained as the balan ce between disk hydrostatic pressure and 
dissociating radiation dHever et alj 120041 ) . Smaller clouds form 
in the outer disk and in some interarm regions, which might have 
a higher SF efficiency and a shorter dispersion time than else- 



Table 1. Radially averaged KS indices (n) and Pearson correla- 
tion coefficients (r). 



SFR 


Molecular 


gas 


Atomic 


gas 


Total j 


jas 


«H 2 






r H, 






r H,„, 


Ha 


1.3 ±0.2 


0.87 


5.6 + 0.9 


0.76 


3.6 + 0.3 


0.92 


FUV 


1.1 ± 0.1 


0.94 


4.1 + 0.8 


0.73 


2.9 ± 0.2 


0.95 


Bol. 


1.1 ± 0.1 


0.94 


4.1 + 0.8 


0.73 


2.9 ± 0.2 


0.95 


mean 


1.2 ± 0.1 


0.92 


4.6 + 0.8 


0.74 


3.1 +0.2 


0.94 



where. GMCs in M33 comprise less than 30% of the molecular 
mass, while in our Galaxy GMCs contain 80% of the molecular 
gas, and the fraction of molecular gas mass in GMCs decreases 
with radius. 



4. Radial Kennicutt-Schmidt law 

We first discuss the radial variations of the KS law by analyzing 
azimuthally averaged values of the SFR and gas surface density 
in bins of 0.24 kpc. To trace the SFR, we consider each of the 
three diagnostics derived in the previous section: extinction cor- 
rected Ha (Eq. [3), extinction corrected FUV (Eq. [4), and bolo- 
metric (Eq.|5]l luminosities. In Table[T]we display the KS indices 
for all three SFR diagnostics, with regard to each gas component 
and to the total gas. The KS indices are rather similar if one con- 
siders different SFR tracers; only slightly higher KS indices are 
found when the Her emission is used to trace the SFR. We shall 
discuss this effect in Sect. [6] 

The KS indices (slopes) are steeper for the total gas surface 
densities (Arecibo plus FCRAO) («h 1o1 ~ 3.1). Since the atomic 
gas, the dominant contributor to the total gas surface density, 
has a shallow radial distribution, we find a very high value of the 
KS index when we correlate the SFR with the atomic gas sur- 
face density alone. A tighter correlation with a shallower slope 
is found between the azimuthally averaged SFR and the molec- 
ular gas surface densi ty with = 1.2 ± . 1, from the inner 
disk out to 6 kpc. The lRobertson & Kravtsovl d2007l) simulation 
results are in remarkable agreement with the radially averaged 
values of the indices found for the atomic and molecular gas in 
M33. 

The KS law is graphically presented in Fig.[3j for the molec- 
ular, atomic, and total gas as a function of the SFR surface den- 
sities calculated from the three diagnostics described above. The 
gas depletion timescales for a constant SF efficiency are repre- 
sented by dotted lines. One can see that the depletion timescale 
for the molecular gas is rather constant and between 0.5 and 
1 Gyr. On the other hand, the depletion timescales for the atomic 
and total gas vary widely across the disk of M 33. The fastest de- 
pletion timescales (about 0.5 Gyr) occur near the center of the 
galaxy where the SF is more pronounced; conversely, the lowest 
depletion timescales are found in the outer parts of M 33 with 
values reaching roughly 10 Gyr. 

Slightly higher indices than what we find using the FUV and 
bolometric emission have been found bv lHever et alj d2004l) by 
examining radially averaged quantities and a SFR traced by IR 
emission alone (e.g., w h, ~ 1.36). Our KS i ndices are compati- 
ble with those found bv lWong & Blitzl ([2002) if we consider the 
molecular gas alone . On the contrary, they are much higher than 
IWong & B litz (2002) for the atomic and total gas, because then- 
sample of galaxies is dominated by molecular gas, while M33 
has less than 10% of its gas in molecular form. 
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Fig. 3. Radial KS law, elliptically averaged over bins of 0.24 kpc. The SFR is calculated from three different tracers: extinction- 
corrected Ha (left), extinction-corrected FUV (center), and bolometric (right). The KS law taking into account only the molecular 
gas surface density is displayed with plus symbols. The KS law relative to the atomic gas is shown by open triangles. Filled circles 
corresponds to the KS law when the total gas surface density is considered. The dotted lines represent the gas depletion timescales 
(cf. labels in Gyr) considering a constant SF efficiency. 



5. Local Kennicutt-Schmidt law 

In this section, we first examine the KS law locally at a spa- 
tial scale of 180 pc (45 arcsec), which is the resolution of the 
FCRAO CO 7=1-0 dataset, and then at larger spatial scales. 
As in the previous section, we calculate the SFR using the three 
different tracers and evaluate the correlations with the atomic, 
molecular, and total gas surface densities. Extinction corrections 
are applied to the Ha and FUV fluxes at each beam position, 
and the SFR(Ha) and SFR(FUV) are computed as described in 
Sect. 12.41 At all spatial scales, there is much dispersion in the 
Ssfr - £ gas relations, and the dynamical range of the gas surface 
density is rather small. As a consequence, the KS indices depend 
on the fitting method used. We shall give the results from two dif- 
ferent fitting methods, which we shall describe in the remainder 
of this section. 

5.1. Local KS law from 180 to 1440 pc: first fitting method 

To estimate the indices of the KS law at a spatial resolution of 
180 pc, without taking into account outlier positions with very 
high or low SFR in our map, we used a recursive fit. The first 
linear fit includes all the 6400 spatial positions de fined by the 
centra l square (80 x 80 positions) of the CO map dHever et al.l 
2004). We estimate the dispersion, <x, of the data with respect to 
this first fit and remove all points that are lying outside 2<x before 
attempting a second fit. We repeat this action recursively until no 
points lie outside the 2<x dispersion and then quote the slope of 
this final linear fit. The value of the threshold has been set equal 
to 2cr, and changing this threshold could lead to slightly differ- 
ent estimates of the final KS index. Too small a threshold will 
remove many points at each step and will end up by discarding 
all the data. On the other hand, if the threshold is too large, the 
final set of data will remain the same as the initial one, as no 
points will lie outside the boundaries. Ideally we want to select 
the value of the threshold which leaves as many points as pos- 



sible, representing the bulk of the distribution in the final sam- 
ple while discarding the outliers; a threshold value of 2cr seems 
appropriate for this study (see Figs. [4] and [5]l and gives stable 
results. In order to study the effect of the resolution on the KS 
indices, we used the same method at coarser spatial resolutions, 
by averaging the surface densities over adjacent positions: 360, 
540, 720, 900, 1080, 1260, and 1440 pc. In Tabled we show the 
final slopes obtained at each resolution for the three SFR trac- 
ers, with respect to the molecular, atomic, and total gas. Also in 
Table [2] we show the Pearson coefficients obtained for the raw 
samples of points. Clear trends appear: the most significant is 
that the correlation between log Zsfr an d l°g ^gas gets better and 
better when the resolution gets coarser. This was expected due 
to the stochastic effects which take place at small spatial scales. 
Although less systematic, the Pearson coefficients are in general 
higher when considering the molecular gas, with respect to the 
total and atomic gas. This latter shows the weakest correlation, 
independently of the SFR tracer used for the correlation. 

In Figs. |4] (molecular gas) and [5] (total gas) we display the 
final slope obtained (heavy line) with the final samples of po- 
sitions included in the fit, for two cases. Positions discarded 
(shown in grey) during the several steps of the iteration process 
lie outside the 2cr boundaries (depicted by dashed lines). One 
can see that the predominant influence on the KS indices is the 
nature of the gas that is taken into account: the molecular gas 
gives values of «h 2 between 1 and 2, while atomic and total gas 
exhibit higher values: n « 2 - 4. As in the azimuthally averaged 
relation, the Ha SFR tracer displays higher KS indices than the 
FUV and bolometric SFR tracers; this will be discussed in the 
next section. 

Indices get steeper on average as the resolution gets more 
coarse, i.e. as we average quantities over larger areas. For the 
total gas density the KS indices obtained radially are even 
higher than the ones obtained for the coarser local resolution. 
Steeper indices for radially averaged densities are also found by 
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IWong & Blitzl d2002l) and iThilker et all d2007l) . For the atomic 
gas the local KS indices are significantly lower, because the H i 
has a shallow radial falloff, while its distribution in the disk is fil- 
amentary. Hence local variations of the atomic gas are more sig- 
nificant than radial variations. This is not the case for the molec- 
ular gas, which declines radially with a scal elength of 2 kpc, 
similar to that of the Ha and FUV emission dPaper III) . The lo- 
cal KS indices relative to the molecular gas are similar to those 
found for azimuthal averages in radial bins. 

The number of iterations needed to reach the final slope 
and the data dispersion generally decreases as a coarser spa- 
tial resolution is considered (in Table [2] the value of the linear 
Pearson coefficient r with the original distributions of points, 
i.e. before the first iteration, is given). This is because stochas- 
tic effects weaken the relation between SFR and gas densi- 
ties at small spatial scales. In the limit of H n region sizes the 
correlation could disappear altogether, because massive stellar 
winds and supernovae explosions could remove the molecular 
gas and quench subsequent SF Using the combined BIMA- 
FCRAO dataset dRosolowskv et al.ll2.007h at a spatial resolution 
of 30 arcsec (120 pc) around selected H n regions, we find in fact 
an even weaker correlation between Xsfr and Eh 2 than over the 
180 pc spatial scale. 

5.2. Relative errors: second fitting method 

One caveat of the procedure described in the previous paragraph 
is that we have considered all positions in the map to allow 
smoothing and averages over larger areas. This implies that even 
positions where the molecular gas detection was not above the 
noise level have been considered. Thus, the reliability of the re- 
sults on the smallest spatial scale relative to the large-scale ones 
could be compromised because of low signal-to-noise; as the res- 
olution is degraded, the signal-to-noise ratio increases because 
we are averaging over larger areas. 

We now fit the log Xsfr — log Sg as local relation at the lowest 
spatial resolution (180 pc), taking into account the errors in the 
determination of the gas surface density as well as in the pho- 
tometry and extinction correction for the determination of the 
Esfr- We use the GALEX FUV map to determine the photomet- 
ric error In the smoothed GALEX FUV image the pixel-to- 
pixel 1-cr noise varies between 2.9 and 8.7xl0~ 6 counts s _1 , and 
we shall consider its maximum value. The maximum large-scale 
sky 1-cr variation of sky uncertainty is 3xl0" 4 counts s _1 . The 
sky noise is much smaller than the large-scale sky variation and 
can be neglected. To the photometri c errors we add the errors 
for extinction corrections. Following Calzetti (2001), we com- 
pute the extinction in the GALEX FUV band as in Eq. [2] The 
largest source of uncertainty in the formula are not the photo- 
metric errors on TIR and FUV luminosities, but the bolometric 
correction factor (1.68). Since the sampled areas contain a mix 
of young and old stellar populations, it is not clear how much of 
the TIR luminosity is due to the young stellar population power- 
ing the FUV emission. We estimate the uncertainties in Apuv to 
be on average 20%. The individual errors on CO flux estimates 
can be computed from the rms noise in individual spectra. To 
these we add the error on the atomic hydrogen surface density, 
which is on the order of 0.7 M pc~ 2 (as from map noise), if the 
total gas surface density is considered. 



2 It is not possible to determine the errors for the Ho- map since the 
map we have is in emission measure units and it is not possible to re- 
cover the original errors in counts. 
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Table 3. Pearson coefficients (r) and slopes (n) derived from the 
"relative errors" method for correlations between the two quan- 
tities shown in the first column. The SFR per unit area, Zsfr, 
is traced by the FUV emission corrected for extinction. Only 
positions where the CO luminosity is above the 2<x limit have 
been considered at a spatial resolution of 180 pc. The surface 
density of molecular, atomic and total gas are Eh 2 > ~^H t and Zh 10 , 
respectively, p is the volume density of the ISM and Ti60 the dust 
opacity. n\ is the OLS slope as described in the text. 



r n m 



log 


ZsFR 


- l°g£nb 


0.42 


2.22±0.07 


1.46 J 


:0.34 


log 


ZsFR 


- log Zh, 


0.24 


1.62+0.03 


1.11=1 


:0.15 


log 


ZsFR 


- log Z H , 01 


0.43 


2.59±0.05 


1.67d 


:0.25 


log 


ZsFR 


-l0g(Z H2 /prSM) 


0.62 


1.16+0.04 


l.OOn 


:0.01 


log 


ZsFR 


- logPlSM 


0.71 


1.07+0.02 


0.97+ 


0.01 


log 


ZsFR 


- logT 160 


0.81 


1.13+0.02 


1.01 + 


0.02 



We now consider all areas in M 33 for which the CO flux 
measured by the FCRAO survey is above 2cr uncertainty. We 
analyze the possible linear correlation between the log of molec- 
ular gas surface density and the log of SFR per unit surface as 
traced by the FUV emission corrected for extinction. A linear 
least square method with errors both i n the gas surfa ce density 
and in the SFR (subroutine fitexy in iPress et al.lll992l) gives a 
slope of 2.22 + 0.07. The Pearson linear correlation coefficient is 
0.42. The error on the slope is determined considering extinction 
correction error variations: 20% on average +5%. The best-fit re- 
gression slopes and correlation coefficients relating SFR and gas 
content in M33 for the 180 pc regions using this method are 
given in Table [3] For comparison, the slopes obtaine d with the 
ordina ry-least-square (OLS) method advocated by Isob e et aTl 
(1990) are also listed in Table[3] This method does not take into 
account either uncertainties or outliers, so is less sophisticated 
than either of the two methods described previously. Indeed, the 
OLS slopes are shallower than either of the other two estimates; 
in the case of poor correlations with many data points, not con- 
sidering properly errors or outliers can give too much relative 
weight to statistically insignificant data. In Fig. [6] we show the 
log £h m - log Ssfr relation with the relative best fitting linear 
relation (KS index = 2.59 + 0.05). The correlation with total gas 
is of similar quality to that with the molecular gas alone, since 
both Pearson linear correlation coefficients are ~0.43. 

We now consider the correlation between the total gas sur- 
face density and the SFR per unit surface, including points which 
have the CO brightness below the detection threshold. For these 
we shall consider only the H i surface density. The correlation 
coefficient does not change, but the regression is shallower with 
a slope of 1 .7 + 0.2. But caution is called for in interpreting such 
a shallow slope as real or claiming that there are areas in the 
galaxy where the surface density decreases but the SFR does not 
(bi-modal distribution). The non-detection of molecular emis- 
sion generates artificial tails with zero slopes at low surface den- 
sities and flattens the average slope. It is more feasible to con- 
sider positions without a reliable CO detection only when the 
corresponding molecular surface density threshold is well below 
the H i surface density, but there are not many positions in our 
data which satisfy this condition. 
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Fig. 4. Local KS law involving the extinction-corrected Ha luminosity as SFR tracer and molecular gas surface density at four 
spatial resolutions. The black points (inbetween the 2cr boundaries, evidenced by the two dashed lines) are the ones considered for 
the last iteration of the first fitting method. The solid line is the resulting best fit considering these positions. The grey points were 
discarded during the various iterations of the fit. The iKennicut t ( 1998b)'s index is depicted by a dash-dot line (slope of 1.4). The 
dotted lines represent the times of gas depletion (cf. labels in Gyr in the lower right panel) considering a constant SF efficiency. 



Table 2. Local KS law indices (n) obtained using the first fitting method for the three SFR tracers. The number of iterations taken 
before the best-fit slope is stabilized is given in brackets. The Pearson coefficients (r) are the ones obtained with the original 
distributions of points (i.e. before the first iteration). 
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180 


1.44 (43) 


0.21 


1.12(24) 


0.24 


1.12(30) 


0.24 


2.53(14) 


0.19 


1.96 (26) 


0.23 


1.96(23) 


0.23 


2.59(18) 


0.21 


1.91 (19) 


0.26 


1.90(19) 


0.26 


360 


1.54(16) 


0.26 


1.17(25) 


0.29 


1.17 (28) 


0.29 


2.59(12) 


0.24 


2.10(14) 


0.29 


2.11 (15) 


0.29 


2.82 (8) 


0.27 


2.09(17) 


0.32 


2.09 (18) 


0.33 


540 


1.64 (20) 


0.29 


1.17(9) 


0.33 


1.17(11) 


0.33 


2.88 (11) 


0.26 


2.18 (10) 


0.30 


2.24(15) 


0.31 


3.39(16) 


0.32 


2.40(13) 


0.38 


2.41 (12) 


0.38 


720 


1.59 (7) 


0.31 


1.16(11) 


0.32 


1.15(9) 


0.32 


3.16(11) 


0.31 


2.44(10) 


0.36 


2.44(13) 


0.36 


3.38 (10) 


0.28 


2.29 (9) 


0.33 


2.30(12) 


0.33 


900 


1.74(10) 


0.38 


1.28 (6) 


0.40 


1.28 (6) 


0.40 


3.29(10) 


0.32 


2.09 (7) 


0.39 


2.10(8) 


0.39 


3.49 (9) 


0.44 


2.22 (6) 


0.53 


2.24 (7) 


0.53 


1080 


1.81 (15) 


0.45 


1.27(13) 


0.55 


1.22 (8) 


0.55 


3.62(6) 


0.32 


2.80 (7) 


0.38 


2.85 (8) 


0.38 


3.98 (13) 


0.48 


2.27 (4) 


0.56 


2.28 (4) 


0.57 


1260 


1.62 (8) 


0.56 


1.18 (4) 


0.63 


1.17 (5) 


0.64 


2.77 (4) 


0.35 


2.82 (5) 


0.40 


2.83 (7) 


0.40 


4.08(11) 


0.52 


3.05 (11) 


0.60 


2.88 (12) 


0.60 


1440 


1.68 (3) 


0.46 


1.38 (8) 


0.51 


1.39 (6) 


0.51 


2.79 (2) 


0.36 


2.73 (3) 


0.46 


2.74 (3) 


0.46 


3.18 (3) 


0.50 


2.60 (4) 


0.61 


2.61 (4) 


0.61 



6. Different tracers of the SFR and of the likelihood 
of the gas to form stars 

In this last section we analyze whether the SFR correlates better 
with other physical quantities of the ISM, such as the hydro- 
static pressure or gas volume density, rather than with the total 



gas surface density. We also discuss the uncertainties related to 
the use of the Ha as SFR tracer, which at first seems more ap- 
propriate than the FUV or bolometric surface brightness in trac- 
ing the most recent star formation. And in the last paragraph we 
show the tight correlation between the SFR and the dust opti- 
cal depth. We briefly discuss whether the dust optical depth can 
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Fig. 5. Local KS law involving the extinction corrected FUV luminosity as SFR tracer and total gas surface density at four spatial 
resolution. The black points (in between the 2cr boundaries, evidenced by the two dashed lines) are the ones considered for the last 
iteration of the first fitting method. The solid line is the resu lting best fit considering these positions. The grey points were discarded 
during the various iterations of the fit. The Kennicuttl (Il998bl) 's index is depicted by a dash-dot line (slope of 1.4). The dotted lines 
represent the times of gas depletion (cf. labels in Gyr in the lower right panel) considering a constant SF efficiency. 



trace the gas prone to form stars in the M 33 disk better than the 
CO 7=1-0 line intensity. 



6.1. The role of the stellar disk and the dependence of the 
SFR on the gas volume density 

The widely used KS law relates the SFR per unit surface with 
the gas surface density. The SFR depends upon the amount of 
molecular gas able to collapse, fragment and form stars on a 
timescale ?sf- Often the relevant timescale has been referred to 
as the cloud collapse or free-fall timescale, which is inversely 
proportional to the square root of the gas volume density in 
the cloud: ?sf fff oc P H ' 5 - If Ph, scales with some power 
of the gas surface density or is constant, then it is conceiv- 
able that the SFR correlates with the surface density of the 
self-gravitating gas, £h 2 , to some power. However, since the 
free-fall time is mu ch shorter than the molecular cloud forma- 
tion t i mescale (e.g..lWard-Thompson et al.lll994tlMcKee & Tar] 
I2003t lEngargiola et dj2003l lKramhol z et al.l2009t) . the relevant 



time for star formation could be linked to the interstellar medium 
out of which the clouds form. For example iKovama & Ostrikerl 
(2009) have recently shown using numerical simulations of disks 
that the ISM structure plays a primary role in determining the 
actual SFR in galactic disks. To test this hypothesis in M 33, we 
relate the star formation timescale not to the molecular cloud 
free-fall time, but to their formation timescale. If t$p is the time 
needed to form self-gravitating molecular clouds out of the dif- 
fuse interstellar gas, f §F will depend on the dens ity of the ISM 
Pism, i e. f s f oc Pjs^ (IKovama & Ostrikerl I2009). The value of 
Pism in the disk scales with the hydrostatic pressure, which in 
galaxies with a low gaseous cont ent such as M 33 depends on 
the stellar disk surface density (ICorbellil I2003 ). An important 
finding in M 33 is that both the outer SF threshold radius and 
the abundance of molecules can be explained if one includes the 
contr ibution of the stellar disk in the h ydrostatic pressure equa- 
tion (lElmegreenl 1 993bt ICorbellill2003h . 

In an attempt to improve the local correlation between the 
gas and the SFR in M 33, we consider the gravitational compres- 
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Fig. 6. Local KS law involving the extinction corrected FUV lu- 
minosity as SFR tracer and the molecular gas surface density in 
(a), the total gas surface density in (b). The spatial resolution is 
180 pc. The best fitting linear relation leads to a KS index of 
2.22 + 0.07 in (a) and of 2.59+0.05 in (b) shown by the heavy 
lines. 



sion of the gas due to the stellar disk. If the hydrostatic pressure 
P sets the gas vertical scale height h g of the disk, the volume den- 
sity pism can be evaluated as P/Cg or following ICorbellj {2003) 
as: 



Pism - 



IK 



where 



ttG \ C„ C s I 



(7) 



(8) 



and S s , c s and c g are respectively the stellar surface density, the 
velocity dispersion of the stars and the gas perpendicular to the 
disk. If pism determines the SF timescale through the cloud for- 
mation timescale (or through ts if pn 2 scales with Pism), we ex- 
pect the following relation to hold: 



I SFR oc 



'SF 



PlSM 



(9) 



We shall consider a gas velocity dispersion of 6 km s _1 
and a ratio of gas-to-stella r velocity dispersion equal to 0.3 
dCorbelli & Walterbosll2007l) . By looking at the correlation be- 
tween Ssfr and Sh 2 /pjoj^ we ^ n< ^ that this is tighter than that 
between £sfr and the gas surface density alone. The Pearson 
linear correlation coefficient is 0.62, and the slope of the linear 
fit is 1.16 + 0.04 in the log - log plane (see Fig. EJ, close to linear 
as expected by the above equation. It seems clear that the gravity 
of the stellar disk, dominant in M 33, plays a major role in reg- 
ulating the SFR. Variations of the gas surface density alone are 
insufficient to explain the scatter in Esfr and imply a steep KS 



Fig. 7. The left-hand panel shows the linear correlation between 
log Esfr and log £h 2 /Pism ( m un i ts °f Mo/pc~ 2 cm -15 ). The 
Pearson linear correlation coefficient is 0.62, and the slope of the 
linear fit is 1.16 + 0.04. The slope is derived taking into account 
uncertainties on both axes, but they are not shown in the graph 
for clarity. On the right-hand panel we show the linear correla- 
tion between logSsFR and logpisM with the Pearson coefficient 
0.71. 



index as well as a loose correlation when considering only the 
gas surface density. 

Finally, it is noteworthy that the Pearson linear correlation 
coefficient is even higher, r = 0.71, when considering the 
log Ssfr - log Pism relation. The slope of the linear correlation, 
shown in Fig. [7] is 1.07 + 0.02. The tighter relation between 
Esfr and pism (or the hydrostatic pressure since pism x P) is 
well known. It can be interpreted as follows: suppose that ?sf 
is proportional to pj^f because the cloud formation timescale is 
proportional to p~££ and regulates the SF. If the molecular gas 
surface density (or volume density if SF is confined into a layer 
of constant scale height), which is in the form of bound, self- 
gravitating units, is proportional to PfS, the expected relation is 
that we find Zsfr k Pism- This would be the case for a constant 
star formation efficiency. 

The surface density of molecular gas in the form of bound, 
self-gravitating units can be different from £h 2 since 2h 2 in- 
cludes the non negligible contribution of diffuse molecular gas. 
Unfortunately there are no detailed surveys of the molecular gas 
in M 33 available as yet which would be able to constrain the 
fraction of gas in the form of bound, self-gravitating units. Only 
results of GMCs surveys are available, but the molecular mass 
spectrum in M33 is steeper tha n -2 and hence dominate d by 
molecular clouds of small mass dBlitz & Rosolow skv 20051). 

6.2. Is Ha luminosity a good SFR tracer on a local scale? 

We compare the three different SFR tracers in Fig. [8] which 
shows Ssfr(Hq;) plotted against Zsfr derived from the FUV (left 
panel) and bolometric luminosity (right) (see Sect. |2). At high 
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Fig. 8. Ssfr derived from Ha plotted against Ssfr(FUV) (left panel) and SsFR(bol) (right panel). Only the 180 pc regions are shown. 
The dotted line corresponds to equality in the two SFRs. 



SFR densities, Ha and FUV give quite similar values, but Ha 
underpredicts SFR relative to FUV at low SFR. SFR as inferred 
from bolometric luminosity is almost always larger than that 
from Ha, and at low SFR densities the effect is quite strong. 

With the Ha emission as SFR indicator, the correlation of 
Ssfr with Z gas or pism becomes looser, and steeper indices are 
found. This is to be expected when one samples regions, as we 
do, with a total bolometric luminosity below 10 39 erg s _1 , and 
we explain why. The bolometric luminosity of a single cluster or 
of an ensemble of clusters measures the mass of the clusters only 
when the IMF is fully populated. Hence, when Lboi is compara- 
ble to or lower than the luminosity of the most massive star in 
the cluster, the IMF cannot be fully sampled up to its high mass 
end. 

As ICorbelli et alj (12009) pointed out, there are two types of 
models for populating the stellar clusters when the IMF is in- 
completely sampled. Either the IMF is truncated to a limiting 
mass which depends on the cluster mass, or the IMF is incom- 
plete and stochastically sampled, but maintains its original shape 
and completeness up to its high mass limit over the whole galaxy. 
For both cases the Lb i/LH ff ratio rises when the most massive 
stars are lacking because the recombination-line luminosity, Lhq- 
has a steeper dependence on the stellar mass than the bolometric 
luminosity does. For the first model, the truncation model, a de- 
viation from a simple scaling law between log Lboi and log Lh is 
expected when Lboi is on the order of the luminosity of the most 
massive star. For 120 M this is about 7 10 39 erg s _1 . For the sec- 
ond model, the randomly sampled cluster model, the stochastic 
character of the IMF implies that low luminosity clusters can be 
made either by populating the IMF up to a certain stellar mass, 
or by just one single star, or by something in between. The pres- 
ence of outliers makes the average L\, ifLu a deviate less from 
its constant value at high luminosities (where the IMF is fully 
populated). In practice it is only for Lboi < 1.5 10 39 erg s _1 that 
Lboi/Lecf drops dramatically. In this regime Leu- is no longer an 



adequate measure of the cluster mass, and for this model a larger 
scatter is expected in the bolometric luminosity - cluster mass 
relation. 

Figure|9]shows the theoretical Lboi - Lhq. relation, referred to 
as the cluster birthline. The filled square symbols correspond to 
the randomly sampled cluster model, and filled circles indicate 
the predic tion for c lusters modelled with a truncated IMF (see 
above and Pape r Hit) . Cross symbols in Fig. [9] show the data for 
the 180 pc regions in M33. The luminosities of these regions 
are clearly sufficiently low to enter the regime of an incomplete 
IMF. These data cannot be used to judge which of the two cluster 
models applies, because regions 180 pc wide ma y contain more 
than one cluster and not necessarily young (see Corbel lTet al.l 
2009, to address this issue correctly). 

The fast decrease of the Ha luminosity in regions of low 
luminosity, which do not contain massive clusters, implies that 
using a constant Ha luminosity - SFR conversion factor, the 
inferred SFR will be much lower than the effective one. If the 
inferred SFR in low density regions is lower, the index of the 
KS relation will be artificially higher. This explains the steeper 
indices we find when using Ha as the SFR indicator in the local 
study of M 33. We conclude that the lack of massive stars in low 
luminosity regions makes Ha an unreliable tracer of the SFR. 

The birthline is the line where very young stellar clusters lie. 
Aging or leakage of ionizing photons bring the clusters above 
the birthline. Aging decreases Lu a faster than Lboi and moves 
Lboi/Lno. to higher values more or less vertically above the birth- 
line (see lPaper IlJ) . Also leakage of ionizing photons from Hn 
regions raises the Lbd/Lp^, values above the birthline. All our 
sampled regions (M 33 areas at 180 pc resolution) lie above the 
birthline (cross symbols) as they should, implying that we have 
not underestimated their bolometric luminosity. However, most 
of the SF regions we sample are well above the birthline, and 
this suggests that there is a mix of young and aged H n regions 
in each area and/or that there is leakage of ionizing photons. The 
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Fig. 9. The theoretical birthline (Paper III) which predicts a de- 
creasing Ha flux for decreasing luminosity of young stellar 
clusters due to incompleteness of the IMF is shown with filled 
squares for the randomly sampled cluster model. The filled cir- 
cles are the model prediction for clusters with a truncated IMF. 
Cross symbols refer to M 33 data, areas observed at 180 pc res- 
olution. 



aging and leakage processes and the approximate formula used 
to estimate Lboi and extinction correction to Lh^ explains the ob- 
served scatter in the Lboi - Lh» relation. M 33 is known to have 
a high diffuse fraction of Ha emission du e to the leakage pro- 
cess f rom individual star-forming regions dHoopes & Walterbosl 
2000). But also i n the continuum UV radiation M 33 has high 
diffuse fractions (Thil keret alj|2005l) . and if radiation escapes 
from the disk from these diffuse regions or from individual star- 
forming sites, both non-ionizing and ionizing photons will con- 
tribute to this leakage. Hence although leakage can explain a 
possible steepening of the Ssfr (from all tracers) at low gas sur- 
face densities, it does not explain the non-linear scaling between 
L bo i or L FUV and L Ha . 

6.3. The correlation between the SFR and the dust optical 
depth 

We now use the dust opacity at 160 fim, Ti6o derived as in 
Paper II, as an unbiased indicator of the gas co lumn density. 
In M33 the metallicity gradient is v ery shallow dMagrini et alJ 
l2007bt iRosolowskv & Simon 120081) . and hence it is probable 
that the gas-to-dust ratio does not change much going radially 
o utwards. The declining behavior of the dust-to-gas ratio found 
in Pa per III when using the 21 -cm emission line brightness to in- 
fer the H i column density and the CO 7=1-0 emission line 
brightness to infer the H2 column density, can then be due to 
an underestimation of the gas surface density at large galacto- 
centric radii. In particular the CO-to-H2 conversion factor can 
be different from the value assumed in the M 33 FCRAO survey 
(2.8 x 10 20 crrT 2 (K km s -1 ) -1 ), and can vary across the M33 
disk according to local variations of the metal abundance and 




-5 -4.5 -4 -3.5 

lQ g T 160 

Fig. 10. Correlation between the SFR (FUV) surface density and 
the dust opacity at 160 fim. Cross symbols are data for areas in 
M 33 which have a 180 pc wide diameter; only those points are 
shown with the CO luminosity above 2cr detection threshold. 



molecular clouds proper ties. Many papers in the literature (e.g. 
iLerov et aill2007l 12009, and references therein) show that the 
CO-to-H2 conversion factor is effectively higher in low metallic- 
ity galaxies and can vary even in single molecular cloud com- 
plexes according to the cloud self-shielding conditions. As an 
alternative there could be opaque H 1 gas unaccounted for by the 
usual optically th in assumption for t he 21 -cm line emission, as 
found in M 3 1 bv lBraun et ail d2009l) . 

We show in Fig. [10] the correlation between the SFR sur- 
face density and the dust opacity. The Pearson linear correla- 
tion coefficient is much higher (r = 0.81) when considering 
the correlation in the logZsFR - logTi6o plane than when us- 
ing the brightness of the 21 -cm and CO 7=1-0 lines to derive 
logXsFR - l°g^H,„, relation (r = 0.43). As shown in Fig. [TOlthe 
distribution seems bimodal; there is an inflection in the slope 
around logri6o ~ -4. The average slope (given in Table and 
the correlation coefficient do not depend on whether we include 
positions where the CO luminosity is below the 2<x limit. We 
have assumed a 20% error on T160 to determine n, but given the 
high Pearson coefficient, the slope is not much dependent on the 
relative errors. 

Although the correlation log £sfr ~~ l°g T i6() is quite tight, 
some more work is needed to check the assumptions made and 
better address the implications of this result. For example even 
though the average radial variation of the metallicity is small 

( 0.4 dex), metallicity rises in the center dU et al .1 12009). and 

at a given radius the variation in the oxygen abundance can be al- 
most one order of magnitude (~ -0.8 dex) (Magrini et al. 2007a) 
. Right now the relative errors on metal abundances are large, 
but if future measurements confirm this scatter, the assumption 
of a constant dust-to-gas ratio, implicit in the correlation, may 
not be appropriate. Also, the actual estimate of the optical depth 
is based on the intensity of the emission at 160 fim and on the 
dust temperature derived using the 70 and 160 yum emission. If 
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the emission at 70 jum is due to dust heated not only by the in- 
terstellar radiation field but also by radiation from star-forming 
regions, then we would expect a higher colour temperature. This 
would give a larger thermal intensity, and thus we would derive 
a smaller T\eo, anti-correlated with the SFR. High-resolution IR 
observations at longer wavelengths, which will be available in 
the near future with Herschel, will allow a more accurate esti- 
mate of the dust optical depth. The resulting correlation coeffi- 
cient involving log Ti6o and log Zsfr could have a higher value 
than the one derived here. 

7. Summary and conclusions 

We test the KS law in the Local Group spiral galaxy M33, 
in azimuthally averaged areas (240 pc wide) and locally from 
~1.5 kpc spatial resolution down to a resolution of 180 pc. 
Starting from a multiwavelength set of observations, we used 
the Ha, FUV, and bolometric (FUV+TIR) luminosities to esti- 
mate the SFR. For gas surface density we consider the molecular, 
atomic, and total gas phases. We use extinction-corrected SFR 
even though extinction marginally affects the observed proper- 
ties of the optical or ultraviolet emissions in M 33 because of the 
rather low dust content of the galaxy. The most important results 
are summarized below: 

- At every spatial scale we find that the Ha KS indices are al- 
ways higher than the FUV and bolometric ones, and we ex- 
plain this as due to the lack of Ha emission in low luminos- 
ity regions where most of stars form in small clusters with 
an incomplete initial mass function at their high mass end. 
We use the cluster birthline to support this, which implies a 
non-linear relation between the Ha and bolometric luminosi- 
ties. The birthline also shows that most regions, even at the 
highest spatial resolution, contain a mixture of ages which 
on average are smaller than 10 Myr. 

- For azimuthally averaged values, the depletion timescale for 
the molecular gas is radially constant and the KS index is 
1.1+0.1 for bot h FUV and bolo metric SFR tracer, lower than 
that found bv lHever et al.l (120041) using IR SFR tracers alone. 
The depletion time for the molecular gas is relatively con- 
stant, with a value of about 1 Gyr. The correlations with the 
molecular and total gas are tighter than with the H i gas, de- 
spite the fact that most of the gaseous mass of M 33 is in H i 
form. Denser filaments of neutral gas are found where star 
formation takes place, but their H i surface density does not 
correlate with the SFR. The KS radial index for the total gas 
surface density is 2.9+0.2. The depletion time for the total 
gas density widely varies from ~0.5 Gyr in the center of the 
galaxy to about 10 Gyr in the outer parts. 

- Locally the dispersion in the SFR-gas density relation per 
unit surface is high, and results are very sensitive to the sta- 
tistical method used to fit the data. The scatter in the KS re- 
lation increases as the spatial resolution increases. At 180 pc 
resolution, fitting the logEsFR ~ logEH 2 data to a straight 
line using a recursive fitting method which removes outlier 
points, we obtain the same the KS index than in the radial 
average analysis (1.1). There is a slight general tendency for 
the local KS laws to show higher slopes for coarser spatial 
resolutions. As for radial averages, the local KS index for the 
molecular gas is much lower than that for the total gas. 

- For our highest resolution (180 pc) we also derive the best 
fitting straight line to the data in the log - log plane, taking 
into account the errors in the determination of the gas surface 
density as well as the errors in the determination of the SFR 



as traced by FUV emission. Considering only the positions 
where the CO detection is above 2<x noise, a bivariate regres- 
sion gives the KS index for the molecular gas of 2.22 ± 0.07 
with a Pearson linear correlation coefficient of 0.42. We be- 
lieve that given the large scatter in the log Ssfr - log 2e 2 rela- 
tion at the highest spatial resolution examined in this paper, 
this method gives more robust result. This KS index is higher 
than the one we find using the recursive fit and all positions 
in our maps. It is also higher than what the ordinary-least- 
square method gives for the same set of data (no uncertain- 
ties considered). The bivariate regression gives a slope for 
the log Ssfr - log £h m relation only slightly higher than for 
the molecular gas alone: 2.64 + 0.07. 

- Given the rather poor correlations between the gas and the 
SFR per unit surface, we analyze whether the SFR corre- 
lates with other physical quantities. A good correlation is 
found with the hydrostatic pressure, i.e. with the interstellar 
medium volume density considering a constant sound speed. 
This implies that the stellar disk, gravitationally dominant 
with respect to the gaseous disk in M 33, plays a major role 
in driving the SFR. The slope of the correlation is close to 
unity, suggesting that the SFR per unit area (or per unit vol- 
ume if the thickness of the SF disk does not vary) is propor- 
tional to the ISM volume gas density. Since the correlation 
is tighter, its slope is less dependent on the statistical method 
used. 

- There is a good correlation between the dust optical depth 
at 160 jum and the local SFR density. This can be inter- 
preted in terms of a Esfr - £gas relation if the dust opac- 
ity is used as an unbiased indicator of the surface density of 
gas prone to star formation. This might be the case for ex- 
ample if the gas-to-dust ratio were constant in M 33 while 
the CO-to-H2 conversion factor were not, or if there were 
opaque Hi gas unaccounted for by the usual optically thin 
assumption for the 21 -cm line emission. The tightness of the 
logEsFR - logri6o correlation and the shallower slope com- 
pared to the logSsFR - logSg as relation would then imply 
that the KS law in low luminosity galaxies still holds, but 
the surface density of the gas in the process of forming stars 
has a wider dynamical range than shown by the intensity of 
the 21 -cm and CO 7=1-0 lines. Still, some caution is 
called for in interpreting this result because of some of the 
assumptions made in deriving it. A more detailed analysis is 
required together with a more accurate determination of the 
dust optical depth, which will be possible in the near future 
thanks for example to the Herschel satellite data. 
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